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ERGODIC BEHAVIOR OF LOCALLY REGULATED 
BRANCHING POPULATIONS 

By M. Hutzenthaler and A. Wakolbinger 

Johann Wolfgang Goethe- Universitdt Frankfurt 

For a class of processes modeling the evolution of a spatially struc- 
tured population with migration and a logistic local regulation of the 
reproduction dynamics, we show convergence to an upper invariant 
measure from a suitable class of initial distributions. It follows from 
recent work of Alison Etheridge that this upper invariant measure is 
nontrivial for sufficiently large super-criticality in the reproduction. 
For sufficiently small super-criticality, we prove local extinction by 
comparison with a mean field model. This latter result extends also 
to more general local reproduction regulations. 

1. Introduction. In naturally reproducing populations one usually en- 
counters an average number of more than one offspring per individual. How- 
ever, given nonextinction, classical super-critical branching processes grow 
beyond all bounds. This is unrealistic because of bounded resources. 

An efficient counteraction to unbounded population growth is achieved 
by a population-size dependent regulation of the reproduction dynamics. An 
example is the so-called logistic branching process [15] in which, in addition 
to the "natural" births and deaths in a super-critical branching mechanism, 
there are deaths resulting from a competition between any two individuals 
in the population. In Feller's diffusion limit, this leads to a negative drift 
term which is proportional to the squared population size and prevents the 
population size X from escaping to oo. This results in an eventual trapping 
of X in zero. 

An attempt to combat this extinction is to consider infinite populations 
modeled by a spatially extended version of the logistic branching process, 
with subpopulations living in discrete demes arranged in the <i-dimensional 
lattice Z d , and with a (homogeneous) migration between the demes. This 
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leads to the following system of interacting Feller diffusions with logistic 
growth with local population sizes X(i): 



Here, the B(i) are independent standard Wiener processes, m is the transi- 
tion matrix of a random walk on Z, d , and 7 are nonnegative constants 
describing the rates of migration, branching and competition. The constant 
K is called the capacity; it is the ratio of the coefficient of super- criticality 
jK and the competition rate 7. 

Models with competition have been studied by various authors: Mueller 
and Tribe [18] and Horridge and Tribe [11] investigated an SPDE analogue 
of (1), with d = 1 and M 1 instead of Z , and Etheridge [7], motivated by 
the work of Bolker and Pacala [2], investigated system (1) and its measure- 
valued analogue (with Z rf replaced by R rf ). These models also include long 
range competition. We emphasize that our methods make use of the fact 
that the interactions due to competition are solely within the same lattice 
site. 

Using arguments from oriented percolation, Etheridge [7] shows that sys- 
tem (1) and also similar systems with nonlocal competition, when started 
from a spatially homogeneous initial state, do not suffer local extinction [i.e., 
Xt(0) does not converge to zero as t — ► 00], provided the capacity K is large 
enough. On the other hand, it was shown in the same paper by a coupling and 
comparison with subcritical branching (similar as in [18]) that a measure- 
valued analogue of (1) with certain nonlocal competition mechanisms suffers 
local extinction. The question whether lattice-based systems like (1) suffer 
local extinction for suitably small K remained open. In the present paper 
we answer this question in the affirmative for the system (1) (Theorem 1). 
More precisely, we specify a strictly positive constant K such that for all 
capacities K < K system (1) suffers local extinction. The constant K de- 
pends on the rates a, (3 and 7 of migration, branching and competition, but 
is uniform in all dimensions d and migration matrices m. 

The proof of local extinction is achieved by comparing (1) with a mean 
field model associated with (1). It turns out that there is no nontrivial equi- 
librium for (1) if is the only equilibrium state of the mean field model. 
This model is mathematically very tractable; its critical capacity, which 
constitutes the upper bound K, is the unique solution of the scalar equation 
appearing in Theorem 1. 

We construct the maximal process X^°°\ which is the solution of (1) 
entering from infinity at time (Theorem 2). As time tends to infinity, 




= a 




(1) 



+ 1 X t (i)(K-X t (i))dt+ y/2l3X t (i) dB t (t) 
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this process converges monotonically in distribution to the upper invariant 
measure of (1). Exploiting a self-duality of X (Theorem 3), we show in Theo- 
rem 5 that the solution Xt of (1), when started in a shift invariant nontrivial 
initial state, converges in distribution to this distinguished equilibrium. 

Self-duality was used to prove ergodicity by other authors, for example, 
Horridge and Tribe [11] and Athreya and Swart [1]. In the latter paper, 
self-duality was established for the resampling selection model which is the 
solution of (1) where the Feller term ^2f3Xt{i) is replaced by the Fisher- 
Wright term ^2[3Xt(i){l — Xt(i)) and where K < 1. Furthermore, Athreya 
and Swart study a branching coalescing particle model which in Feller's 
diffusion limit leads to the solution of (1). For both models, they prove 
existence of the maximal process and of the upper invariant measure. 

2. Main results. We obtain the local extinction result and the result 
about existence of the maximal process and of the upper invariant measure 
for a more general class of interacting locally regulated diffusions, in which 
both the individual fertility and the branching rate depend on the local 
population size Xt(i). The system of stochastic differential equations we are 
going to consider is 



dX t (i) = a[Y j rn{i,j)X t (j)-X t (M dt 



(2) 

+ h{X t (i)) dt+^2- g(X t (i)) dB t (i), i e G, 

where G is an at most countable Abelian group. Unless stated otherwise, 
we will assume that the migration matrix m is translation invariant, that 
is, m(i,j) = m(0,j — i), that X]jeG m (CM) = 1j an d that m is irreducible, 
that is, 3n:nS n \i,j) > 0. Assumptions on the functions h and g will 
be collected in Assumptions Al, A2 and A3 below. 

An appropriate state space for (1) and (2) is provided by a construction 
going back to Liggett and Spitzer [17]: For given to, let a = {o~i)i^G be 
summable and strictly positive such that 

(3) ^2aim(i,j) < C LS o-j, j£G, 

i 

for some Cls < oo, and define the Liggett-Spitzer space 

(4) E a = |aE€Rf :||£|| ff ■=J2 a i\ x i\ < °°|- 

Notice that every translation invariant measure \i on M> with f XQii{dx) < 
oo is supported by E CT . 

The following assumption guarantees existence and uniqueness of a strong 
Eq- valued solution of system (2). 
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Assumption Al. The functions h : M>o — > M and g : M>o — > K> are lo- 
cally Lipschitz continuous in [0, oo) and satisfy h(0) = g(0) = 0. In addition, 
the function h is upward Lipschitz continuous, that is, 

(5) sgn(x -y)(h(x)-h(y))<C\x-y\ 

for some constant C. Furthermore, g is strictly positive on (0, oo) and sat- 
isfies the growth condition 



(o) hmsup <oo. 

x— >oo X 

Proposition 2.1. Assume Assumption Al. Then, for any i£E ff , i/ie 
system (2) /ias a unique strong solution X = (Xt) starting in x and having 
almost surely continuous paths in the normed space E CT . 

This will be proved in Section 3. Together with Assumption Al, the fol- 
lowing condition on the drift function h will be essential in Theorem 1. 

Assumption A2. The function fr:R> — >-M is concave. Moreover, for 
some xq > 0, h is negative on [xo,oo), and satisfies 

r+oo j 

(7) / — ——dx<oo. 



■>-■<) 



-h(x) 



Note that (7) together with the upward Lipschitz continuity required in 
Assumption Al implies 

(8) lim h(x) = — oo. 

x — >oo 

For the interacting Feller diffusions with logistic growth (1), the functions h 
and g are of the form 

(9) h{x) = jx(K — x), g(x)=(5x. 

In this case, Assumptions Al and A2 are clearly satisfied. 

Theorem 1. Assume Assumptions Al and A2. Denote by X the solu- 
tion of equation (2) for an arbitrarily prescribed initial distribution on E CT . 
If 

, . f°° h(y) ( fy-ax + h(x) , \ , 

(10) / -^" ex P / ' dx)dy<0, 



g{y) \h g(x) 

then the process suffers local extinction, that is, 
(11) (X t )=><So ast^oo. 

Here, denotes the zero configuration. 
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In the logistic Feller case (9), condition (10) simplifies to 

(12) J exp\Kjy - —y 2 j ■ aexp(-ay)dy <1; 

see the proof of Corollary 6 at the end of Section 5. 

A first glance at system (1) might tempt one to believe that, even for small 
capacities K (and a fixed), a suitably mobile migration m in the dynam- 
ics (1) could prevent the system from suffering local extinction. However, 
Theorem 1 and condition (12) reveal that this is not the case. 

The proof of Theorem 1 will be given in Section 6. Its main idea is a 
comparison with a mean field model corresponding to (2), given by the 
solution V of 

(13) dV t = a(EV t - V t ) dt + h(V t ) dt + ^2-g(V t )dB t . 

For this comparison, we will suppose that the initial condition Xq of (2) has 
a distribution ft which is associated, that is, 

(14) J fi-f 2 dfi> J hdflj J 2 dji 

for all bounded, coordinate- wise nondecreasing f±, f 2 :E CT — ► R. 

For the solution X of (2), we will show that, for all t > 0, the marginal 
distributions of Xt are bounded by the distribution of Vt in the <j cv - order 
(where "icv" stands for "increasing, concave"; see [21] for this and related 
notions). More precisely, in Section 6 we will prove the following. 

Proposition 2.2. Assume Assumption Al and concavity of h. Let X 
be a solution of (2) whose initial distribution p, is associated. Assume that 
the Xq(i), i £ G, are identically distributed and have finite expectation. Let 
Vt = (Vt(i))ieG be a system of processes coupled through the initial state 
Vo(i) = Xq(i), i £ G, but following independent mean field dynamics, that 
is, every Vt(i) solves equation (13) with Brownian motion B{i), where the 
B{i), i G G, are independent. Then 

(15) Wf(X t ) < Wf{V t ) 

for all bounded, coordinate-wise nondecreasing and concave functions f : E CT — > R 
depending on only finitely many coordinates. 

The following theorem, whose proof will be given in Section 4, provides 
for the existence of a distinguished equilibrium state of (2), called the upper 
invariant measure. For proving Theorem 2, we will exploit the following 
assumption, which relaxes Assumption A2. Condition (16) ensures that the 
drift is "sufficiently negative" for large values of X t {i) so that the process 
"comes down from oo." 
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Assumption A3. There exists a function h > h such that, for some 
xq > 0, h is negative and concave on [xq,oo) and satisfies 

r+oo \ 

(16) / — s dx<oo. 

Jxq —h(x) 

To prepare for Theorem 2, we need a bit of notation. If Hi,^2 are probabil- 
ity measures on a partially ordered set S, then we say that /ii is stochastically 
smaller than or equal to fj,2, and we write [X\ < fj,2, if there exists a random 
pair (Yi, I2) with marginal laws (Yi) = [/,{, i = 1, 2 and Yl < I2. We say that 
a sequence of probability measures /U n increases stochastically to a proba- 
bility measure ^oo, denoted by pn f ^oo, if there exists a random sequence 
(Yi) which a.s. increases to Y^ and has marginal distributions (Yi) = Mi) 
i = 1, 2, . . . , 00. 

Theorem 2. Assume Assumptions Al and A3. There exists an E^- 
valued process (-X^ )t>o wifft the following properties: 

(a) For eac/i e > 0, (X^°°')t> e is a solution of (2) starting at time t = e. 

(b) Tae /irs£ moment of \\ X^ 00 ^ \\ a is finite for every t > 0. 

(b) Lei x( n ) = (x^^jgc, n = 1, 2, . . . , &e an increasing sequence in E CT sac/i 
iaai, /or aH jgG, 

(17) | 00 as n — > 00. 

If X^ is the solution of (2) starting in x^ £ E CT ai izme zero, i/ien 

(18) C{x[ n) )] £(X t (oo) ) as n 1 00 (t>0). 

(d) There exists an equilibrium distribution v ( called the upper invariant 
measure,) for the dynamics (2) such that 

(19) £(X t (oo) )|P as t too. 

(e) ^4ny E, a -valued solution X of (2) satisfies 

(20) £(A 4 ) < £(X t (oo) ) (f>0). 

/n particular, any equilibrium v is stochastically smaller than or equal to v . 

(f) 5oi/i f/ie upper invariant measure v and C(Xt^) are translation in- 
variant and associated. 

In the following two theorems, we exploit the specific form of the dy- 
namics (1) of the interacting Feller diffusions with logistic growth. As it 
turns out, the solution of equation (1) has a property of self-duality which 
is helpful for the investigation of convergence to equilibria. 
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For the formulation of the self-duality result, write m) for the transpose 
of the matrix m, choose a satisfying (3) with rrv instead of m, and recall 
that E CT t denotes the corresponding Liggett-Spitzer space. 

Theorem 3. Assume f3 > 0. Let X and be solutions of (1) with 
migration kernels m and rrv . Then we have the following self- duality: 

(21) E*exp(-I<X t ,y)) = E^exp (-^(x, 

for all x G K a , y £ E CT t ■ 

A similar (though nonself-) duality for interacting Feller diffusions (or 
super- random walks), that is, (1) with 7 = 0, is given by 

(22) ESexp(-<Xt,y)) = exp(-<x,i%>), 
where v = (yt(i)) solves the initial value problem 

(23) ^v t (i) = J2m(i,j)(v t (j) -v t (i)) -v t (i) 2 , i£G,v = y;, 

see, for example, Chapter 4 of [5]. 

The proof of Theorem 3 is contained in Section 7. The main advantage 
of the self-duality (21) is that instead of starting in a configuration with 
infinite total mass, we can analyze the evolution of the process started with 
finite total mass. For example, choose y = \5o and x with x(i) = const. Then 
the self-duality tells us that it makes no difference whether we study the law 
of X t (0) started in x, or that of the total mass \x}\ :=X)i-^t(*) with X^ 
started in X5o, A > 0. This leads to the following corollary (see Lemma 7.1 
together with Theorem 1): 

Corollary 4. Assume f3,j > 0. Let the parameters a,(3,j,K be such 
that inequality (12) holds. Then the solution X of (1) started from an initial 
state of finite total mass [i.e., Y^i-^o(i) < °°] hits in finite time a.s. 

In the more general situation of (2), self-duality is not available for proving 
the analogue of Corollary 4. We conjecture that also for interacting locally 
regulated diffusions (2), the local extinction (11) implies global extinction 
of finite mass processes as stated in Corollary 4, but we do not have a proof. 

Theorem 3 will be the principal tool for proving convergence to the up- 
per invariant measure specified in Theorem 2. This convergence will be the 
subject of Theorem 5 below. On an intuitive level, the reason for this con- 
vergence is as follows: There are two forces working against each other, 
super-critical branching and individual competition. The third ingredient 
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is migration, which is important for spreading out newly produced mass. 
Super-critical branching increases mass, whereas competition among the in- 
dividuals decreases it. If a (local) population size is large, then competition 
is stronger, whereas, as long as a local population size is small, the compe- 
tition is negligible in comparison to the mass producing branching. Thus, 
there should be some attracting equilibrium state in which the two forces 
balance each other. This is the upper invariant measure. 

Theorem 5. Assume /?, 7 > 0. Let X be a solution of (1) and suppose 
that C{Xq) > fj,, where fi is a measure on E CT which is translation invariant 
and does not charge the zero configuration 0. Then 

(24) £(Xt)=>P as i-> 00, 
where v is the upper invariant measure. 

From this it is clear that the only extremal translation invariant equilib- 
rium distributions are 5o and v. They coincide in case of local extinction and 
differ in case of survival. Section 8 will be devoted to the proof of Theorem 5. 

3. Preliminaries. We turn first to the proof of existence and uniqueness 
of (2) under Assumption Al. Let us define 

(25) bi(x) := ai^2m(i,j)xj - Xij + h(xi), xGE CT , 

\jeG ) 

where a = {*Ji)ieG satisfies (3). Denote z + := z V 0. By inequality (3) and 
Assumption Al, there exists a finite constant C\ such that 

(26) ^^l x ,„ K >o(6 l (x)-6 i (y))<C 1 ||(x.-y.) + || (T Vx,ySE, 

for every subset MCG. From inequality (26), we will obtain monotonicity 
in the initial configuration. This monotonicity is a crucial property which we 
will exploit several times. First, we prove boundedness of second moments. 

Lemma 3.1. Suppose that h and g satisfy Assumption Al. Let (Xt) 
be any weak solution of equation (2) with E||Xo||^ < 00 whose paths are 
continuous in E CT . Then there exists a constant C < 00 such that, for each 
T>0, 

(27) supE||A t ||2 < (1 +E||A ||2) e CT < 00. 

t<T 
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Proof. Let Gk be finite subsets of G which monotonically exhaust G 
as k — > oo. Denote H^Ho-^ := J2i£G k a i\ x i\- Applying Ito's formula, we obtain 

d\\X t \\l k = 2\\X t \\ atk J2 (Ti{bi{X t )dt + V2g(X t (i))dB t (i)) 

ieG k 

(28) 

+ 2 y £o*g(X t (i))dt. 



i&G 



k 



Let n E N. The continuous function g is bounded on the interval [0,ra/o"j] for 
every i S G k - Thus, the stochastic integrals on the right-hand side of (28) 
are L2-martingales when stopped at time r n := inft>o{||Ai|| CT > n}. By path 
continuity, we have T n — > oo as n — > oo almost surely. Taking expectations, 
inequality (26) with y = implies 



E\\ X tAT n \\l,k 

(29) 

<E|LY ||L+2E 



J^ Tn (ci\\X 8 \\l+J2<$9(Xs(i))^ ds. 



r n \\l)ds 



By the growth condition (6), we know that g(x) < C2(l + x 2 ) for some con- 
stant C2 < oo. Letting k — > 00 and using monotone convergence, we obtain 

(30) E\\X tATn f a <E\\X \\l + C 3 f\l + E\\X. 

Jo 

for some constant C3 < 00. Applying Gronwall's inequality to the function 
1 1— > 1 + E||Xi AT - n \\ 2 , we arrive at 

(31) E||X Mr X <(l + E||Xo||^)e C3t -l. 
Letting n — > 00, Fatou's lemma completes the proof. □ 

In the proof of Proposition 2. 1 , we need a stronger uniformity than Lemma 3.1 
provides. 

Lemma 3.2. Assume Assumption Al. Let (Xt) be any weak solution of 
equation (2) satisfying condition (27). Then for each T > 0, there exists a 
constant Ct < 00 such that 

(32) Esup \\X t \\a < C T (1 + E||A || CT + E\\X \\l) < 00. 

t<T 

Proof. Recall the definition of G k and || • from the proof of Lemma 3.1. 
Multiplying by dj and summing over i £ Gk in (2), we obtain, for t <T, 

\\Xt\\ a ,k — ||A"o|| -,fc 

(33) t t 

= / ^bi{X s )ds+ / o-iV2g(X s (i))dB s {i). 
Jo i£G k Jo ieG k 
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The estimate (26) implies that J2ieG k a i^i{^s) < Cill-^slU- Thus, denoting 
the rightmost term in (33) by M/% we obtain 

(34) sup||X w || CT)fc < ||X || CT + / Ci sup \\X r \\ a ds + sup |M*|. 

u<t JO r<s u<T 

The process (M t fc ) is an L2-martingale since, by the assumption g(x) < 
C{\ + x 2 ) and condition (27), the integrands y/2g{X s {i)) in (33) are square 
integrable, and the second moment E|M^| 2 = Jq 2 J2ieG k a i^9{Xs{i)) ds is 
bounded by Cy(l +E||Xo||^.) for some constant Ct- Thus, using the esti- 
mate z < 1 + z 2 , we conclude from Doob's ^-inequality that 

(35) Esup|M^| < 1 + E|M|| 2 < 1 + C T {1 + E||X |£)- 

u<T 

Therefore, taking expectations in (34) and applying monotone convergence, 
we obtain 

Esup ||X U || CT < E||Xo|| CT + C\ [ Esup ||A" r || CT (is 

u<t JO r<s 

(36) 

+ 1 + Ct(1 + E||X || 2 ) 
for all t <T. Now the assertion follows from Gronwall's inequality. □ 

The following monotone coupling lemma will be an important tool. 

Lemma 3.3. Let hi, hi and g satisfy Assumption Al, and let B = (B(i))i^G 
be a system of independent Brownian motions defined on some filtered prob- 
ability space. For i = 1,2, assume that X L is defined on the same probability 
space, satisfies equation (2) with Brownian motions B(i), drift function h L 
and initial configuration i l 6E ff , and has continuous paths in E CT . Then 

(37) h\ < /12 together with x} < x 2 implies X\ < X 2 Vt > a.s. 

Proof. The first part of the proof follows that of Theorem IV. 3. 2 in 
[12] .Let 1 > ai > ■ ■ ■ > a n > ■ ■ ■ > be defined by 

rl \ ra\ \ ran-i \ 

(38) / —du = l, — du = 2,..., —du = n, 

Notice that a n —> as n — > oo. For every n = 1,2,..., define a continuous 
function ijj n (u) with support in (a n ,a n _i) such that 

(39) 0<Vn(n)< — and / ip n {u)du = l. 
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Furthermore, define 

px ry 

(40) (j) n {x):=t x> Q dy ip n {u)du, x£R. 

Jo Jo 

These functions satisfy 4> Tl E C 2 (R), |^(x)| < 1, 4>n{x) = 1 x> Qil) n {x), 4> n {x) < 
x + and 4> n (x)^x + as n —> oo. Fix i G G and let Tk := inf{i > 0:X^(i) V 
X?(i) > k}. Write Af := X t x (i) - X?(i) and let be as in equation (25) 
with h replaced by h L , i = 1,2. By Ito's formula, 

rtATfc 



(41) 



^(A l s )[V2 5 (X s 1 (,)) - V2 5 (X a 2 (i))]d£ s (i) 

+ / <^(A;)[&^i)-fe 2 )]^ 

<(A*)[V2 5 (Xj(z))- V2 5 (X 2 (z))] 2 cfe. 



As n — ► oo, the left-hand side converges to (Aj ATfc ) + — (A ) + in L\ by domi- 
nated convergence and Lemma 3.1. In the rest of the proof, C±, C2, . . . will be 
suitably chosen finite constants. By Assumption Al, there exists a constant 
C\ such that g{x) < C\(l + x 2 ). Thus, Lemma 3.1 implies Eg(X^(i)) < 00 
and we have, by dominated convergence, 

tAT k 



(42) E/ (t A%>0 -4>' n (Al)) 2 (V2g(Xl( i ))-V2g(X 2 m 2 ds n ^0. 



Hence, the first (stochastic) integral on the right-hand side converges in 
L2 to the same expression with 4>' n (x) replaced by l x >o- For the second 
integral, notice that b\ is globally Lipschitz continuous on {x : X{ < k}. Thus, 
for s < Tfe, |6|(Xg)| is bounded by C2||A^|| CT , which has finite expectation by 
Lemma 3.1, and we obtain, by dominated convergence, 

(43) / |l Ai>0 -^(A«)|-|6j(X. 1 )-M)|(fa , ^0. 
Jo 

Finally, we consider the third integral on the right-hand side of equation (41). 
The local Lipschitz continuity of g implies that yfg is globally l/2-H61der 
continuous on the interval [0, k] . Therefore, the last integral in (41) is bounded 
by 

(44) / .C 3 |A*j<^^0 asrwoo. 

Putting these calculations together, equation (41) implies 

(Aj)+-(A )+= / t l A|>0 [V2 5 (X s 1 (i)) - V2g(X 2 s (i))} dB s (i) 
Jo 
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(45) 



t^ibKxD-bfix^ds 



for all t < T/% almost surely. By path continuity, we have r\. — > oo almost surely 
as k —> oo and thus, equation (45) holds for all t > 0. The stochastic integral 
on the right-hand side is an /^-martingale because of g(x) < C\(l + x 2 ) and 
Lemma 3.1. Taking expectations, we arrive at 



E 



5>(A|)+ - $>(A* ) 



i&G 



i&G 



(46) 



<C 4 f E^^AD+ds. 

Jo ieG 

In the last step, we used b\ < b\ and inequality (26). By Gronwall's inequal- 
ity, we obtain 



(47) 



E||(A|)- 



< 



(A J ) + | 



,c A t 



ieG. 



For later use, we note that this inequality implies 



(48) 



E||A||| CT < WA^e 



c 4 t 



i€G, 



if 6i = bo- For this, notice that \x} 



The 



right-hand side of inequality (47) is zero by the assumption x 1 < x 2 , which 
finishes the proof of the monotonicity result for fixed t > 0. Finally, X\ < X 2 



follows for all t 6 
surely. □ 



and then by continuity of paths for all t > almost 



Proof of Proposition 2.1. Let B = (Bi) ie a be a system of indepen- 
dent Brownian motions, and fix an initial condition i£E ff . We will prove 
existence of a solution of (2) similarly as in [9], where the system (2) is stud- 
ied in the case h = 0. To this end, for finite ACG and i,j G G, we define 
m A (i,j) := m(i,j)ljj g A an d consider the finite-dimensional system 

dX^i) = « E m\i,j)X t A (j) dt - aX A (i) dt 



+ h(Xf(i)) dt + V2- g(X^(i)) dB t (i), is A. 

Under Assumption Al, equation (49) has a unique solution X A starting in 
(xj)jgA. We extend X A to an infinite sequence (still denoted by the same 
symbol) by putting X^(i) := for i & G \ A. Following the arguments in 
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the proof of Theorem 1 in [9], one can show that there exists a process 
X = (Xt(i)) arising as the monotone limit 

(50) X^(i)]X t (i) as A t G. 

To show that X has a.s. continuous paths in E CT , we first note that, for each 
finite A C G, the process X A , being a finite-dimensional diffusion, has a.s. 
continuous paths and therefore satisfies 

(51) limPf sup \\X A -X A \\ a >e) =0 

S ~*° \\t-s\<6,s,t<T J 

for all e > and T > 0. 

For all finite A C G, the process X A satisfies the assumptions of Lemma 
3.1, with m(i,j) in (2) replaced by m A (i,j). Consequently, X A also satis- 
fies (27), where the constant C can be chosen uniformly in A. Therefore, 
by the monotone convergence (50), X satisfies (27) and, due to Lemma 3.2, 
also (32). 

Next, we set out to show that, for all e > and T > 0, 

(52) limPf supWXt-XfWa >s) =0. 

A TG \t<T J 

For this purpose, let Gk and llxllo-^ be as in the proof of Lemma 3.1. 
From (50) together with the a.s. component- wise continuity of X and Dini's 
theorem, we conclude that, for all T > and fc£N, 

(53) suppf* ~X^\\^ k ^0 a.s. as A | 67. 

t<T 

By (32) and dominated convergence, we therefore have 

(54) Esup||X t - X^Wvk ^0 a.s. as A | 67. 

t<T 

For every finite ACG and fcGN, we estimate 
Esup\\X t -X A \\ a 

t<T 

(55) 

<Esup||X t -A: t A || (Tifc + 2Esup V OiX t (i). 

The rightmost term in (55) does not depend on A and converges to 0, again 
because of (32) and dominated convergence. Together with (54) this implies 
that the left-hand side of (55) converges to zero, and proves (52). 
For e, S and T > 0, we have the estimate 



P( sup \\X t -X s \\ a >3e 

^\t-s\<S,s,t<T 



(56) 

\t-s\<8,s,t<T / \t<T 



<P( sup \\Xf-X A \\ a >£) +2P( sup||X t -X t A || CT >£ 
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Because of (51) and (52) the left-hand side of (56) converges to as 5 — > 0. 
This implies almost sure pathwise continuity. 

For uniqueness, we proceed as follows. In the situation of Lemma 3.3, 
choose hi = hi and x 1 = x 2 . Then pathwise uniqueness follows by applying 
Lemma 3.3 twice. Uniqueness in law and strong existence follow then from a 
Yamada-Watanabe type argument (see [22], Theorem 2.2). For the existence 
of a strong solution, it remains to show that the dependence of the unique 
solution on the initial configuration is measurable. This follows from the 
monotonicity result of Lemma 3.3. □ 

Lemma 3.4. Let h and g satisfy Assumption Al. The strong solution 
Xt of system (2) is monotonically continuous in its initial configuration in 
the following sense: Let x^ n \x£ E CT , be the starting points of X^ and Xt, 
such that 

(57) x {n) T (i)x as n | oo. 
Then 

(58) A t (n) t (i)X t Vt > as n f oo a.s. 

PROOF. In equation (48), let h 1 = h 2 := h, X\ := X t and X? := x[ n) . 

Letting n — > oo, this implies Li-convergence of Xt — X\ for fixed time t > 0. 
The monotonicity result of Lemma 3.3 finishes the proof. □ 

4. The upper invariant measure. Proof of Theorem 2. 

Proof of Theorem 2. To fix notation, let us write C-{X t ) for the 
distribution of X t [the solution of (2)] starting from an element x € E CT . For 
N G N, we define the element N G E CT by N(i) = N , i G G. Let xf be the 
process started from N_. By Lemma 3.3, the sequence Xf- is nondecreasing 
in iV for all t > 0; let us write X^ 00 ^ for its a.s. limit. 

Now let (x^) be a sequence as in Theorem 2(c). For all n G N, we conclude 
from Lemma 3.4 that 

(59) C{X[°° ] ) \ C^(X t ) >£^ (n)AN (X t ) / C^ n \x t ). 

Again by Lemma 3.4, we obtain, for all iV G N, 

(60) C^ (n)AN (X t ) / C K (X t ). 

n— >oo 

Thus, by a diagonal argument, there is a subsequence xj nN ^ of x^ such that 

(61) C^ {nN)AN (X t ) / £(X 4 {oo) ). 

N^oo 
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Together with inequality (59) and monotonicity (Lemma 3.3), this results 
in 

(62) C^ N \x t ) / £(X t (oo) ). 

As (x^) is an increasing sequence, (62) is equivalent to (18). 

The next step shows that the limit is finite almost surely. Let h>hbe the 
function given by Assumption A3. Notice that h may be replaced by h + C 
for every constant C > 0. Furthermore, h is bounded above. Thus, we may 
assume that h>h\s concave. By Ito's formula, Lemma 3.1 and translation 
invar iance, 

(63) j t EX fW = m*r(i)) < EM^fW) < Hvxf(i)). 

For the last step, we applied Jensen's inequality. Therefore, the expectation 
is bounded above by the deterministic function y(t,x) satisfying 

(64) ^y(t,x) = h(y(t,x)), y(0,x)=x. 

The concave function h(x) converges to — oo as x — > oo. Choose xq such that 
h is strictly negative for all x > x$. Then for all x > xq and t > 0, we have 
< y(t,x) < x. From (64), we obtain by separation of variables that the 
solution satisfies 

rx 1 

t = - dz 

Jy(t,x) h(z) 

( 65 ) ! ,00 I 



< / — dz I / — dz asx^oo. 

Jy(t,x) —h(z) Jlim^^oo y(t,x) —h(z) 

For the monotone convergence, notice that y(t,x) is nondecreasing in x and 
that all integrals are finite by inequality (16). Hence, if lim^^oo y(t, x) was 
infinite for t > 0, then we would face the contradiction < t < 0. Therefore, 
we arrive at 



E\\X^\\ a = ^a^ lim EXf(i) 



(66) 



W-»oo 



< Oi hm y(t, x) < oo, t > 0. 

L — ' x— too 



From Lemma 3.4, it is then clear that, for all e > 0, the solution of (2) which 
starts at time t = e from Xe°°^ is the a.s. monotone limit (as — ► oo) of the 
solutions of (2) starting from Xs~ at time e, or, equivalently, starting from 

N at time 0. At the beginning of the proof we defined x[°°^ as this limit; 
hence, we have so far proved parts (a), (b) and (c) of Theorem 2. 
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A similar argument as in (59) proves that the process with initial measure 
fx is dominated by the maximal process, which is part (e). 
To prove part (d), fix < s < t. By part (e), 

(67) £(X^)>£ x t-s(X r ). 
Using this with r = s, we get the inequality 

(68) £(a1°°)) > £ x(t ~° (x s ) = £p^ (oo) ), 

where the last equality follows from the Markov property. We conclude from 
this monotonicity that C(X^ 00 ^) j v for some probability measure v on E CT , 
which by continuity in the initial configuration (Lemma 3.4) is an equilib- 
rium distribution of the dynamics (2). 

Next, we show that the upper invariant measure is translation invariant 
and is associated. Both properties are preserved under weak limits. Further- 
more, we will argue that these properties are preserved under the dynamics. 
The constant configuration X^f = N is both translation invariant and as- 
sociated. Hence, both Xf- and have these properties for all t > 0. 
Therefore, the claim follows. 

The translation invariance of the migration kernel implies that the dy- 
namics (2) preserve translation invariance. To prove the preservation of as- 
sociated measures, we will argue in a similar way as in [4] where the analogue 
of (2) with h = and [0, 1] instead of R> was treated. We first consider 
the approximation scheme (X A , A) with finite AcG, used to prove the ex- 
istence part of Proposition 2.1. For fixed A, Theorem 1.1 in [10], together 
with a uniform approximation of h and g on compact intervals by smooth 
and bounded functions and with inf- E >ofl , fe(x) > 0, shows that, for an 
associated initial distribution £(Xq), the projections of £{X^) to M> are 
associated. Since £(X^) approximates £(X t ) as A j G, the claim follows. 
□ 

5. The mean field model. In this section we study the dynamics 

(69) dV t = a(EV t - V t ) dt + h{V t ) dt + V2g(V t )dB t . 

It can be shown (but will not be required for the subsequent proofs) that (69) 
arises as the limit of a sequence of processes following the dynamics (2), 
where G is replaced by a finite set G n of cardinality n and rrS- n \i,j) = 
1/n for i,j E G n . This type of limit is known as mean field or Vlasov- 
McKean limit; we will therefore address (69) briefly as mean field model. 
Intuitively, a uniform migration which spreads out mass as far as possible 
should be good for survival, and conversely, extinction of V governed by 
(69) should imply extinction of X governed by (2). With this motivation in 
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mind, we investigate in this section conditions on h and g under which the 
dynamics (69) admit a nontrivial equilibrium distribution. 
To this end, we consider the following: 

Lemma 5.1. Suppose that Assumption Al holds and that 

(70) 3y >0:h\ [0iyo] >0 and j£ % 0jOo) < 0. 

There is no nontrivial invariant measure for the dynamics (69) if and only 
if 

ta , f°° h(y) ( fy-ax + h(x) , \ . ^ n 

(71) / -^-exp / v ' dx )dy < 0. 



\Jyo a{x) 

If condition (71) is not satisfied, then there is exactly one nontrivial in- 
variant measure. Assume additionally that lim^o //° "g^^ dx exists in 
(—00,00]. Then condition (71) is equivalent to 

f°° ay ( fy -ax + h(x) , \ , 
(72) / -r^exp / — ^drU/<l. 

Proof. Let > and consider the process given by 



(73) cttf = a{6 - V t 9 ) dt + h(V t e ) dt + V2g(V°)dB t . 

By standard theory (e.g., pages 220f and 241 in [14]), the equilibrium dis- 
tribution of (73) is 

(74) T e (dy) = ^-e,J f y ^- X ^ h ^ dx)dy=:C e Hy)dy, 

g{y) \Jyo 9{x) J 

where Cq G (0, 00) is the normalizing constant. Indeed, existence of an equi- 
librium of (73) is clear since the drift in zero is positive in zero and becomes 
sufficiently negative near 00; formally, this follows from the finiteness of the 
integral Jq°° &(y) dy, which can be checked easily. 

Obviously, (69) admits a nontrivial equilibrium if and only if f yTo(dy) = 
6 has a positive solution. Hence, all we need to do is to characterize the 
situations where 

(75) $9>0:f(e):=aJ^±r e (dy) = 0. 

We eliminate one occurrence of 6 on the left-hand side of (75) by an inte- 
gration by parts: 

Jo g{y) \J yo 9{x) J \J yo g(x) J 
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(76) = lim 

e-+0 



. v a (0 - x) + h{x) , 
exp / — dx 

yo 9{x) 



+ / ^ fexp (r^dx)) 9 exp ( ~ aX : frfc> dx ) dy. 



g(y)\ \Jy g( x ) J J \Jy g{x) 

We now analyze the two boundary terms on the right-hand side of (76). In 
the following calculations, Cj are finite constants. Recall that h is nonpositive 
for large arguments. Furthermore, in Assumption Al we assumed g{x) < Cx 2 
for some constant C and all x > yo > 0. With this, the expression coming 
from the boundary value 1/e tends to zero as e — > 0: 



(77) 



/ f 1 / 6 a(0-x) + h(x) , 
0<ex P (/ dx 

\jyo gy- 1 ) 

\Jy ve Cx A 



For the other boundary term, we recall that h is nonnegative for small 
arguments and estimate 

/ fv° a(6 -x) + h(x) , 

< exp - / -i — ^ dx 

V Je g{x) 

/ rvoA(e/2) q/2 
(78) <C 2 exp-/ ^4-dx 



g(x) 



V 2 7 + 



dx 

By assumption, (7 is locally Lipschitz in zero and thus, g(x) < C4X in a 
neighborhood of zero. Together with > 0, this implies that all boundary 
terms vanish. Notice that the expression coming from the boundary value e 
does not need to be zero in case 6 = 0. 

At this point we have seen that / can be rewritten as 

( 79) m =rMUjf"" te))^jf"^^ dx ) is 

Jo g{y)\ \Jyo g{x) J J \J yo g{x) J 

for 6 > 0. We will show that / is strictly decreasing and continuous in 9 > 0. 
For this, consider the function 

(80) (op (jT^* 

for fixed y > 0. If y < yo, then h(y) > and the integral is negative. If 
y > yo, then h{y) < and the integral is positive. In both situations, the 
function in (80) is nonincreasing. Furthermore, there is an interval [2/1,2/2] 
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with uq < yi < y2 where h(y) < and where the function in (80) is strictly 
decreasing and converging to — oo. The integral over [0,yo] on the right-hand 
side of (79) is continuous and nonincreasing in 9 > 0, and bounded in 9 > 1. 
This follows from dominated convergence and the fact that the integral over 
[0,e] on the right-hand side of (79) is bounded above by 



(81) 



/ "TT ex P / n~r V ) y SU P h\ x ) ex P ( / r~\ dx 

Jo g{y) \Je 2g{x) J x < yo \J yo g[x) 



v a6 , 

- exp( I Mx~) dx 



£ 



£ 

■ C < oo 





for all 6 > 6 > 0, where e > is such that \ax — h(x)\ < a9/2 for all x < e. 
By monotone convergence, the integral over [yo,oo) on the right-hand side 
of (79) is continuous and strictly decreasing in 6 > 0, and decreases to — oo. 
Thus, the function / is continuous and strictly decreasing in 9 > with 
/(oo) = — oo. Hence, condition (75) is satisfied if and only if lim^o /(#) < 
0. Note that, by strict monotonicity of /, there is at most one nontrivial 
invariant measure. 

For the limit 9 — > in equation (79), we use monotone convergence (for 
the Jq° part) and dominated convergence (for the part). Thus, we have 

(82) hm/ W = f^expf f ZSE+M^. 

e^o w Jo g(y) \J yo 9&) J 

Therefore, ]img^of(9) < is equivalent to condition (71). 

Now, additionally assume that lim e _+o JJ -^r~^ dx exists in (— oo, oo]. 
Then reversing the calculation in (76) with 9 = 0, we arrive at 

iim/w= r^expf r z^+m dx ] dy 



+ o g(y) \J yo g(x, 

(83) 

/ . fy° -ax + h(x) , 
— exp — hm / — ax 



s- 



'oJ e g(x) 



If the limit on the right-hand side is oo, then lim^o /($) > anci a nontrivial 
invariant measure exists. The assertion is true in this case because the left- 
hand side of (72) is oo. Otherwise, the limit on the right-hand side of (83) is 
finite. Then multiply the equation with exp(J^° dx) and merge the 

two integrals Jq° and J^ Q into one integral. Hence, we see that (71) and (72) 
are equivalent. □ 



We now specialize this result to the logistic Feller case, where condi- 
tion (72) can be simplified. 
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Corollary 6. Consider the mean field model (69) with h(x) = jx(K — 
x) and g{x) = fix. Assume a, 7,/3 > and let K > be uniquely determined 
by 

(84) J exp(Kjy-—y 2 yaexp(-ay)dy = l. 

There is no nontrivial invariant measure for (69) if and only ifO<K< K . 

Proof. First of all, convince yourself that Assumptions Al and A2 
hold. Thus, Lemma 5.1 applies if K > 0. After an integration and a change 
of variables (y — » (3y), condition (72) takes the form (12). The left-hand 
side in (12) is strictly increasing in K, tends to oo, is continuous in K by 
monotone convergence and is smaller than one for K = 0. Hence, K exists 
and is unique. By monotonicity, condition (12) holds if and only if K < K. 
□ 

For example, in the case a = 7 = (3 = 1 formula (84) gives the numerical 
value K = 0.6973 

The following extinction result for the mean field dynamics is a fairly 
direct consequence of Lemma 5.1. 

Lemma 5.2. Consider the mean field model given by (69). Suppose that 
Assumptions Al, A3 and condition (70) hold. Then inequality (71) implies 
local extinction: 

(85) £(y t )=>5 (as t-> 00) 
for any initial law. 

Proof. Paralleling the arguments in Section 4, one infers the exis- 
tence of the maximal process V^ 00 ^ for the dynamics (69), which obeys 

Ciy^ 00 ^) > CiVt). Again, this maximal process converges to an invariant 
measure. However, by Lemma 5.1 and condition (71), the trivial measure <5o 
is the only invariant measure. This implies the assertion. □ 

6. Comparison with the mean field model. Proof of Theorem 1. The 

main idea for the proof of Theorem 1 is the assertion that the interacting 
locally regulated diffusions are dominated by the mean field model. The 
intuition behind this is that a uniform spread of mass reduces competition 
and therefore is good for survival, and that the mean field model arises as a 
limit of uniform migration models (see Section 5). 

We proceed in two steps to prove Theorem 1. First, we establish a com- 
parison between the system of interacting locally regulated diffusions (2) 
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and the mean field model (69) which implies that it is more likely for the 
latter to survive. Then we exploit the fact (proved in Section 5) that for 
some parameter configurations not even the mean field model survives. 

The proof of the comparison result will first treat the case where the 
functions h and g satisfy the following assumptions. 

Assumption A4. The set J is a closed finite interval of the form [0,c], 
< c < oo. The functions h : I — ► R and yj~g : I — ► R are twice continuously 
differentiable on / and satisfy h(0) = g(0) = g{c) = > h(c). Furthermore, g 
is strictly positive on (0, c). 



The proof of Proposition 2.2 is based on the following lemma. 



Lemma 6.1. Let h and g satisfy Assumption A4. Suppose that h is 
concave and that the set A is finite and nonempty. Then the semigroup of 
the solution of equation (49) preserves the function cone 

(86) P = {/ 6 C? l(R |„) :£,>„ <0 VU}, 

where C^ 1 (R> ) denotes the space of all bounded C 2 - functions /:R> — > R 
with bounded first partial derivatives. 

Proof. This lemma is an addendum to Proposition 17 in [3]. There, the 
preservation of F was proved for h = and matrices m with J2jeA m (hj) = 1 
for all i € A. This proof also works for more general matrices m which only 
satisfy J2jeA m (hj) < 1 for i G A. To extend the argument to the case h^O, 
let y(t,x) be the solution of 

(87) ^y(t,x) = h(y(t,x)), y(0,x)=xel. 

This defines a deterministic Markov process whose semigroup is given by 
Stf(x) := f(y(t,x)). Similarly as in [3], we only need to establish that this 
semigroup preserves F if h is twice continuously differentiable. A little calcu- 
lation shows that it is enough to prove that y(t,x) is increasing and concave 
in x. To show concavity, notice that differentiating equation (87) results in 



d d 2 ( d \ 2 d 2 

For fixed x, write (88) as z' t = at + bf zt with zq = 0. The solution for this is 



Zf = exp^y b s ds^j J exp^— J b r dr^ja s ds. 
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Since h(x) is concave, at is negative, implying the claimed concavity. A 
similar, even simpler argument shows monotonicity. □ 



Proof of Proposition 2.2. We make use of the approximation scheme 
X A defined in the proof of Proposition 2.1; recall that X A is the solution of 
(49). Since A f A j X t , it suffices to show the inequality (15) with Xt replaced 
by Xf~, and / depending only on the coordinates X{ with i & A. 

Furthermore, we assume for the rest of the proof that h and g satisfy 
Assumption A4. The general case follows then by approximating h and g 
pointwise by functions hk and gk satisfying Assumption A4. See Lemma 19 
of [3] for the details. 

In addition, we may assume that / G F; otherwise approximate / by 
functions in F and use dominated convergence. Denote by St the strongly 
continuous semigroup of A A defined on C(I A ). When applied to (p E F, the 
generator of X A takes the form (see, e.g., Theorem 7.3.3 of [19]) 

By Proposition 1.1.5(c) of [8], we know that 

(9i) j t s t f = gs t f. 

Let Vt = (Vt{i))i£\ be a system of processes coupled through the initial 
state Vq(i) = Xq(i), i G A, but following independent mean field dynamics: 



(90) Qtp{x) 



E 

ieA 



a 



<p(x). 



df(V t ) 
(92) 



a E j-f(7t)(&V t (i) - Vtii) + h(V t (i))) dt 



ieA 



ieA 



dxi 



where B(i)i^\ are independent Brownian motions. Write fit := C(Vt); for 
brevity, we suppress in this notation the dependence on A. By equation (92), 
the evolution of fit is given by 



(93) 



ieA 



E' 



{WVt{i)-Vt{i))[^f){V tj 



d d^ 1 



Integration by parts yields 



(94) 



Jo (^Vs^St- s fds = [fi s S t _ s f] t -J^ Hs-^Stsfds. 
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In view of (91), this reads as 

(95) fkf - fiStf = (j^P's- fisQ^jStsfds. 

We will show that the integrand is nonnegative. Prom Lemma 6.1, we know 
that ip := St-sf (for < s < t fixed) is an element of F. By equations (90) 
and (93), 



(96) 



d _ _ r \ 

= a^E^Fi(f)E^ 

ieA 



-a m{i,j)W 
> a^E' r F(i)E' 2 

ieA 







<p)(Yt) 



-a £ m{i,j)WV t (j)W 

i,j£A 

Note that under the assumptions on p, we have C(V(i)) = C(V(0)) for all 
i G A. The right-hand side of (96) is nonnegative because of X^eA m (hj) ^ 1- 
To see the inequality in (96), notice that is bounded and component- 

wise increasing by Lemma 6.1. The claimed inequality thus follows from the 
fact that C(Vt) is associated, which we now prove. Independent real-valued 
random variables are associated (see page 78 of [16]), and p is associated by 
assumption. Hence, 

mf(V t )g(V t )} = E*[E[/(F t ) 5 (F t )|F ]] 

>E^[E[/(F t )|F ]E[ 5 (F t )|F ]] 

(97) > E*[E[/(F t )|Fo]]E^[E[ 5 (Ft)|F ]] 

= EP[f(V t )}W[g(V t )}, 
showing that C(Vt) is associated. □ 

Proof of Theorem 1. Let y := max{y > 0:h(y) = 0}. Assume for 

the moment that yo > 0. The measure p := C(x[°°^) is associated, shift 
invariant and its first moment is finite by Theorem 2. Let (Vt) be the solution 

of (69) with initial distribution fi := C(x[°°\o)). Theorem 2(e), (a) and 
Proposition 2.2 imply 

-\X t+1 (i) > J3 e - Ax t+1 (i) _ -£p t 



(98) 



-AX t (oo) (i) > E M e -A 1 
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It follows from Lemma 5.2 that, under the stated assumptions, Ee~ 1 — ► 1 
for all A > as t — > oo. This proves the assertion for the case yo > 0. 

If yo = 0, then h(x) := (h(x) — h(l)) A satisfies h < h < because h is 
concave. Let X be the solution of (2) with h replaced by h and with the 
same family of Brownian motions. By the previous step, X suffers local 
extinction. Lemma 3.3 implies X < X , which completes the proof. □ 



7. Self-duality. Proof of Theorem 3. 



Proof of Theorem 3. In the rest of the paper we exploit the specific 
form of the dynamics (1) for the interacting Feller diffusions with logistic 
growth. Theorem 3 states that the process is "dual to itself" via 



(99) 



~(X t ,y) 



>,4: 



We will prove this for the solution X A of (49). By (51), we know that the 
process X A monotonically approximates X. Hence, the assertion follows by 
dominated convergence. 

For the rest of the proof, we consider X A . We write X instead of X A 



and x,y instead of x A ,y A . The duality function is H(x,y) 



Recall the definition of C^( 



exp(-^(x,y)). 



from Section 6. Define the linear operator 



C(R A ) by 



Q x f{x) 



(100) 



E 








J dxi 



+ 'yx i (K ■ 



d 

Xi)^—f + (3x t 

OXi 



2-t 



By Ito's formula, the process (Xt)t is a solution of the martingale problem 
)). In order to apply Theorem 4.4.11 of [8] (with the choice 



61 



for {g x ,c 

a, (5 = 0), we will show that 



>o 



(101) 



G X H(;y)(x) 



G Xt H(x,-)(y) 



Vx,yE 



»A 
l >0- 



We prove equation (101) by considering the different parts of (100) sepa- 
rately. Since H is a function p((x,y)) of the scalar product, it is easy to see 
that the migration terms of both sides are equal. To establish equation (101), 
it remains to show that 



lXi{K 



d d 2 

X i) — H ( X ' y) + x i^-2 H {X, y) 



(102) 



dxi 



dxi 

d d 2 
Vi)-^-H(x,y) + f3yi—^H(x,y) 
oyi dyf 
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for all i £ A. Observe that this equation is symmetric in x and y. Consider 
the left-hand side of equation (102) divided by H(x,y): 

jXi{K -Xi) ■ (-Zvi) + ^ Xi ' y H^) 

(103) 

1 2 K 7 2 2 7 2 2 

= -n-Xiyi + —xfyi + -jXiVi ■ 

The right-hand side of (103) is symmetric in x and y and therefore, by 
interchanging the roles of x and y, equals also the right-hand side of equa- 
tion (102) divided by H(x,y). 

Theorem 4.4.11 of [8] is applicable if we prove that 

(104) sup \g x H(x s ,x})\ 

s,t<T 

is integrable for all T < oo, where X and X^ are independent. It is not hard 
to see that 

(105) \g x H(x,y)\ <C(\x\\y\ + \x\ + \y\) Vx,yGM| 

for a finite constant C. For this, use that zexp(-z) is bounded in z > 0. 
Integrability of (104) therefore follows from the independence of X and X^ 
and from Lemma 3.2. □ 

Let us write Ai c (Z d ) for the set of configurations in ]R>q with finite sup- 
port. As a consequence of the self-duality, we prove the following character- 
ization of the upper invariant measure in terms of the finite mass process. 

Lemma 7.1. Assume /?, 7 > 0. The upper invariant measure v of (1) is 
uniquely determined by 

J exp(-^{x,\)jp(dx) = P A (3t>0 such thatxj = 0), |A| < 00, 

(106) 

where X' is the solution of (1) with the transpose migration matrix rrv . 

Proof. Fix a configuration A £ M. c {'L d ) and consider the process (X^) 
started in the constant configuration n(i) = n. This process converges to the 
maximal process as n — > 00. Therefore, the self-duality implies, for t > 0, 

Eexp(-2(X<~ ) ,A>)= n lim o E e xp(-i{4" ) ,A)) 



E A lim expf 1 ■ 

n— »oo 



(107) 

= P A (A/ = 0) 

= P x (3s<t:Xl = 0). 
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For the second equality, we used monotone convergence. Letting t — > oo, 
the assertion follows from Theorem 2(d). For general A with |A| < oo, use 
monotone convergence. □ 

8. Convergence to the upper invariant measure. Proof of Theorem 5. 



Proof of Theorem 5. Let \x be a translation invariant distribution on 
7L d with /x(0) =0. For analyzing the long-term behavior of the interacting 
Feller diffusion with logistic growth started in fx, we can assume without 
loss of generality that \i has finite first moment and satisfies (jl{xq = 0) = 0. 
Otherwise, we let the system run for a little time e > 0, obtaining 



(108) 



lim C^Xt) = lim C c ^ Xe \x t ). 



t— >oo 



t— >oo 



A comparison with the maximal process [see Theorem 2(e), (b)] yields 
E^X e (0) < oo. Furthermore, after a fixed positive time e > 0, every com- 
ponent is strictly positive almost surely (see Lemma 8.2). 

Let X and X^ be solutions of (1) with migration matrix m and its trans- 
pose m' , respectively. In Lemma 8.1 we will show that the total mass hits 
zero in finite time or tends to infinity. Hence, we get by self-duality (Theo- 
rem 3) 

7, 



(109) 



E»exp(-j(X t ,\) 
H(dx) ( E A 



x 



exp 



We treat the two terms on the right-hand side separately and begin with 
the first term. Apply Holder's inequality to the integral with respect to /x. 
For this, let 1/pi = Xt(i)/\Xt\ if this is positive. Thus, we obtain 



E ; 



l \xt\~*oo 



(110) 



E ; 



1, 



fi(dx) exp 

Yl [ I /i(dx)exp 



■^X-IXM 



fi(dx)exp[ -—x \X. 







Xt(i)/\Xt\ 



as t 



oo. 



The equality is a consequence of the translation invariance of \x. The last 
expression tends to zero because of dominated convergence and the assump- 
tion /i(xq = 0) = 0. As to the second term on the right-hand side of (109), 
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dominated convergence gives 



li(dx)E x [l 



3s:X]=0> 



(111) J v(dx)B x \l 3s:Xt= ^ W (-j(x,x}) 

as t —> oo. Using Lemma 7.1, we arrive at 

Hm E^exp^-^(X t ,A)^) = P A (3t >0 s.t. xj = 0) = J exp 

Starting in C(Xq) > fx, the process C(X t ) is bounded below by C^(X t ) 
(Lemma 3.3) and is bounded above by £(X t ) [Theorem 2(e)], which both 
converge to v. This concludes the proof of Theorem 5. □ 







(x, A) I v(dx). 



We have to append the following: 



Lemma 8.1. Assume f3 > 0. Let J 6e a solution of (1) starting inx £E ff 
u«t/i finite total mass \x\ < oo. Then with probability 1, either, 

• t/iere is a t>0 such that X s = /or all s>t or 

• >oo ast^oo. 

Proof. The intuition behind this is the following. The process always 
has a positive probability of hitting the lower trap. Whenever the total mass 
stays bounded, the process will seize its chance. 

This is made precise in Theorem 2 of [13]. In order to apply this result, we 
only need to verify that there always is the risk of extinction in the following 
sense: 

(112) Vy: inf P x (3t : X t = 0) > 0. 

\m\<y 

Let 1^ be a solution of (2) with G = Z d , h(x) = jKx and g(x) = fix. By 
Lemma (3.3), X and Y may be coupled such that Xt is bounded above 
by Yt almost surely. Furthermore, |Yj| is equal in distribution to Feller's 
branching diffusion Ft with super-criticality jK started in \x\. The extinc- 
tion probability of Ft is strictly positive; see, for example, Appendix 6.2 of 
[6]. Therefore, condition (112) follows from 

(113) P x (3t >0:X t = 0)> Pl x l (3t > : F t = 0) > P y (3t > : F t = 0) > 
for every x with \x\ <y. □ 

Lemma 8.2. Suppose that h and g satisfy Assumption Al. Let X be 
a solution of (2). If its initial law \i is translation invariant and does not 
charge the zero configuration 0, then, for every fixed time to > 0, 



(114) 



X to (i)>0 VieGP^-a.s. 
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Proof. Assume that h < 0. Otherwise, compare Xt with the process 
defined with h A instead of h. 

Let = — a.9 + h{9). For e > 0, define the solution of 



(115) dY t e,i = aedt + h(Y t s,i ) dt + y/2g(Y t £>i ) dB\, Y^ > 0, 

on the same probability space as X by using the same system of Brownian 
motions. This system satisfies P°(Y^' 1 > 0) = 1 for all t > 0. Otherwise, 
continuity in the initial value would imply that there is a t > and a 9q such 
that pOtYt' 1 = 0) > for all 9 < 0q. Integrating this with the equilibrium 
distribution Y £ [see equation (74), it exists because of h < 0] yields r e (0) > 0, 
which is false. Thus, we have 

y£m(i,j),i q 

(116) 5 

VeG(0,l)nQ V5G(0,l)nQ Vi,j s.t. m(i,j)>0 a.s. 

Denote the event {X t (j) > eVt G [to — 5, *o]} by ^4 e>< 5. On A e _s we compare X 
with the solution of (115): 



(117) 



X t {i) = X to _ s (i) + [ a^2 m (i,k)X s (k)ds 

+ f* h(X s (i))ds+ f V2g(X s (i))dB s 

Jtn-8 Jtn-S 



ho— 8 Jto—8 

> / am(i,j)eds + / h(X s (i))ds 
Jt -S Jto-8 



+ /* V2g(X B (i))dB 8 (i), 

Jtn-8 



»t 
lto-8 

for all t € [to — 5, to] . By standard comparison results (e.g., Theorem (V.43.1) 
in [20] and a stopping argument), this implies Xf (i) > yT 1 on A £ s a.s. 
By path continuity, A £t $ approximates {Xt (j) > 0} as 5, e — > 0. It follows 
that on X to (j) > we have X to (i) > for all i such that m(i,j) > a.s. 
With the migration kernel being irreducible, every site can be reached from 
j. By induction, we conclude that every component of X to is positive a.s. 
given X to (j) > 0. 

Starting in a nontrivial translation invariant measure, the system a.s. 
never hits 0. Therefore, there is a location j with Xt (j) > a.s. This proves 
the lemma. □ 
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